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1. Introduction 

Nonlinear Shrodinger equation (NLSE) and its vector modifications have been studied 
in various contexts for more than 30 years. It has found applications in such areas 
as plasma physics, hydrodynamics, magnetism and nonlinear fibre optics. The NLSE, 
the so-called modified NLSE, and some of their vector counterparts (e.g. Manakov 
equations) are integrable and, as such, possess stable soliton solutions that exist due 
to the precise balance between the effects of nonlinearity and dispersion. Since solitons 
preserve their shape and internal properties during the propagation they have been 
widely used in fibre optics as information carriers PJ E] . 

In most applications, however, integrable systems represent only the leading 
approximation of the real physical model. Quite often such models can be analysed 
using the concept of a nearly integrable system. In such systems the deviations from 
the initial integrable model are considered to be small and the underlying dynamics can 
be inferred via methods of perturbation theory [Hj. The important and distinct class of 
perturbations constitute random perturbations of different origin. In particular, in fibre 
optics, stochastic perturbations emulate the effects of the amplifier spontaneous emission 
as well as various fibre inhomogeneities (HEUHj. The influence of random perturbations 
on solitons has been studied in different integrable systems including Korteveg-de-Vries 
(KdV) equation [HEHi, NLSE jHH El EE2 ECU EE! (see the latter Ref. for additional 
references), modified NLSE [TI)J El EI, and Manakov system [3J El . Inasmuch as 
solitons in integrable systems are robust against small perturbations, the weak noise 
leads to random walk of soliton parameters (which is called jitter) as well as generation 
of dispersive radiation. In the sequel we will only be concerned with the additive white 
Gaussian noise (AWGN) model, that represents the simplest physically justified example 
of random perturbation. For such systems most celebrated result for soliton jitters 
is the so-called Gordon-Haus effect: the random walk of the soliton position due to 
the effect of weak AWGN [10J. The main tool used for studying the statistics of the 
soliton parameters for the aforementioned systems is the adiabatic soliton perturbation 
theory (see e.g. review paper [3] and references therein). Since the perturbation has 
stochastic nature the resulting system of equations for soliton parameters represents a 
Langevin system of stochastic equations with the multiplicative noise. For this Langevin 
system, in the vast majority of works on the subject the noise was assumed to be 
additive, and the equations were linearised. The assumption of the additive noise in 
the linearised Langevin equations can give rise only to the Gaussian probability dencity 
functions (PDF). However the exact statistics of the NLSE soliton, for example, in the 
presence of AWGN are not Gaussian as was shown in papers [TTJE1E1E!- Gaussian 
approximation can only be used for approximate estimates of the integral characteristics 
of the PDF, like e.g. variances of soliton jitters. 

In the present paper we will consider an integrable Manakov system driven by 
vector AWGN, a system which, in some respect, is a vector generalisation of the NLSE. 
Soliton of the Manakov equations possesses two additional parameters as compared to 
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NLSE soliton. Doktorov and Kuten fTI)] demonstrated in the Gaussain approximation 
that timing jitter of the Manakov soliton can be managed by adjusting these two extra 
parameters. In the current paper we will show that not only timing jitter but all 
the remaining jitters can be managed in the similar fashion provided that the vector 
AWGN has anisotropic statistical properties. We will also obtain (for the first time 
to our knowledge) marginal PDFs for the parameters of the Manakov soliton, that 
describe correctly the non-Gaussian statistics of rare fluctuations. We will derive the 
exact statistics of the soliton parameters using the combination of the adiabatic soliton 
perturbation theory and the Fokker-Planck equation (FPE) approach, proceeding in the 
same way as in Refs. fSl H3J for the NLSE soliton. 

The paper is organized as follows: in Section El we derive a system of Langevin 
equations with multiplicative noise using well established technique of the adiabatic 
perturbation theory for a single Manakov soliton. In section El the Fokker-Planck 
equation attached to the corresponding Langevin system is derived. In the same section 
we find the exact solution for the PDF of the amplitude- frequency jitter. Then we use 
a WKB approximation for the FPE to obtain a marginal PDF for the soliton amplitude 
and polarization angle. Eventually, in Section H] the Gaussian limit is considered - 
we calculate variances for all the soliton parameters. In Conclusion we summarise the 
results obtained. 



2. Stochastic adiabatic perturbation theory for Manakov soliton 

Consider a vector soliton propagating along a random medium, which is described by 
the perturbed Manakov equations 

i* z + *« + 2(***)* = n, (1) 

where the components of vector ^ = (u(z, t), v(z, t)) T represent two polarization states 
of a propagating vector field (e.g. electrical field in a birefringent optical fibre), whereas 
the r.h.s. of Eq.([TJ accounts for the random noise (e.g. ASE) and its components have 
the following correlation properties: 

(n a (t, z)) = (n a (t, z)n/3(t', z')) = , (2) 
(n a (t, z)n* tf, z')) = D Qf3 5(z - z')5(t -t'), a, (3 = 1, 2. (3) 

Here matrix D is symmetric and positive definite. 

In the absence of any perturbation the Manakov equations belong to the class of 
those which are integrable by the inverse scattering transform method. The general form 
of the one-soliton solution of Manakov equations ^ is given by following expressions 
(see e.g. 0): 

u(z, t) = 2r](z) cos f3{z)sec\i{2ri(z) [t - T(z)}} exp{ia(z) - 2i [t - T(z)\ £(z) + i(f)(z)}, (4) 
v(z, t) = 2r){z) sin /5(2)sech{2^(^) [t - T{z)}} exp{ia(z) - 2i [t - T(z)} £(z) - i<f>(z)}, (5) 
and dependencies of the six soliton parameters on the propagation distance z are: 

T){z) = T] , (6a) 
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£{z) = 6>, (66) 

= 4 (^o +^o) 2; + «o, (6c) 

T(z) =-4£ z + T , (Qd) 

ip(z) = <po, (6e) 

= A), (6/) 

where the constants with zero subscripts are the values of soliton parameters at z — 0. 

Because the noise intensity is always small, the noise term, n(t, can be considered 
as a perturbation. In the presence of an arbitrary perturbation, R = (Ri,R 2 ) T , the 
soliton parameters become slowly varying functions of propagation distance z and their 
evolution during propagation should be determined via the equations of the adiabatic 
soliton perturbation theory 

dry 

dz 

dz 
da 
cL7 
dT 
cU 

^ = R e y dt g ;(t,z)R(t,z), (7e) 

^ =R e y dtg;{t,z)R{t,z), (If) 

where vector functions gj can be expressed through derivatives of the single-soliton 
ansatz *o with respect to appropriate soliton parameters. These functions can be 
written in the compact form as follows: 

19*o 1 9* 

g " ~ 2~dc7 ' ft - ~- &» + ^-fff ' 

csc(2/3)9* 19*o 

- cos(2/3)g v + 



Re y dtg^(t,2)R(t,2), (7a) 

Re J dtg^(t,z)R(t,z), (7 b) 

4(n 2 + e) + Re j dtg*(t,z)R(t,z), (7c) 

-4£ + Re / dtg^(t,z)R(t,z), (7d) 



gr 



4?7 9/3 ' v ^ JOV 2n d£ 2 dn 

1 9*o cot(2/3) csc(2/?)9* 



Ar] 9£ ' 2?7 &r ' At] dip ' 

The domain of integration in (fTaj)-j7fl) should be from — oo to oo. The set of equations 
(J7aJ) -(|7j|) with R(t, z) = n(t, z) describes the fluctuations of soliton parameters under 
the action of AWGN provided that the interaction between the soliton and background 
radiation is negligible. The perturbation theory is not applicable in the limit 77 — >- and 
also in the limit when the polarization angle (3 approaches values 7m/2 with an integer 
n. 
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3. Marginal PDFs for soliton parameters 

In this section we proceed to derive the FPE for all six soliton parameters and obtain 
the exact statistics of the amplitude-frequency jitter. We also use a WKB approach to 
derive marginal PDF for soliton amplitude and polarization angle. For simplicity we 
consider the case where the components of perturbation vector n are two independent 
noises with the same intensity D: 

(n a {t,z)np{t',z'))=D5 a p5(z-z')5(t-t'), a,p = l,2. (8) 

In that case the FPE reads as (see Appendix): 



dP 

dz 



dr/ 2 

df3 2 
d 2 



Dt] 
' D 



16 7] 



P 



P 



+ 



+ 



da 1 



Dn 
~12 
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+ 



d 2 

dT 2 



Dn 2 
192 r] 3 



-P 



+ 



_l_ + vrV + 3a 
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dipda 



£> cot (2/3) csc(2/3) 



lQi] 



P 



d 2 



dTda 



Drr 2 ^ 
96n 3 



+ 



P 



8P_ 

dip 2 

cot 2 {2(3) 



Dcsc(2pf 



16 7] 
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+ 
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Dcot(2/9) 
877 



P 



d_ 

da 



(9) 



The boundary conditions in angular variables a, tp and (3 are those of periodicity. In £, 
r) and T the PDF must decrease rapidly enough to provide normalization. At r] = we 
assume that the rj component of the probability density current vanishes (see below) 
which means that no solitons are created of annihilated. 



3.1. Marginal PDF for amplitude-frequency jitter 

One can obtain in a closed form an autonomous FPE for the PDF of the amplitude- 
frequency jitter P(r],^\ z). Thin can be done by integrating out the redundant degrees 
of freedom in (JHJ) by virtue of boundary conditions. Then the FPE for the PDF of the 
amplitude- frequency jitter is 



dP 

dz 



d_ 

dn 



3D 



p 



+ 



drf 



Dn 



p 



+ 



Dr) 
~12 



P 



(10) 



Since, without loss of generality, the initial value of £ may be put equal to zero the latter 
equation must be solved under the initial condition 

P(v^\0) = 5(OS(v-Vo)- 
The //-component of the probability current vector is zero in the plane 77 



:iv, 



0: 





d 


'Dn 


) 






drj 






T)=0 



0. 



(12) 



The solution of this equation can be obtained by virtue of Fourier transform with respect 
to £ and Laplace transform with respect to n. The appearing equation is then solved 
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by the method of characteristics. Omitting the details, the solution is found to be 



P(v^\z) = — I Pk(r]\z')e-^dk . 

Z7T 



where the Fourier component is written as 
P k {r]\z ) = — ^r-—^ — - exp[-a fc (?7 + rj ) coth(a k z )\ I 2 1 



(13) 



(14) 



7/0 sinh(a/ c z / ) ' \smh(akz') 

Here I 2 (. • •) is the second-order modified Bessel function, a k = k/y/3 and z' = Dz/A. 
Integrating Eq.(|13|) over £ we obtain the marginal PDF P(i]\ z): 



i j) r 

P(V\z) = -— exp <^ - 



V + Vo 



z' 



It is easy to find that the asymptotic of the PDF for the amplitude jitter is 



P( V \z) 



n 



3/4 



exp 



2V^V /4 

The higher-order momenta of amplitude and frequency are: 



oo . 



(v n ) 

,, d 2n 

Q k 2n 



|sech 3 (a fc z') exp[— a fc r/ tanh(a/ c 2; / )] } 



n 



1,2, etc. 



(15) 

(16) 

(17) 
(18) 



k=0 



The odd momenta of frequency vanish. In equation ()17)) for the amplitude momenta 
L^(. . .) are the generalized Laguerre functions. In particular, one can obtain the growth 
of the average amplitude with the propagation distance: 



3D 

W =Vo + —z. 

3.2. Marginal PDF for polarization angle and amplitude. WKB-solution 
Marginal PDF P(r),0\z) satisfies the following FPE: 



(19) 



OP 

dz 



d_ 



d_ 

dp 



D cot {2(3) 
8^ 



P 



+ 



drf 



Drj 

T 



p 



+ 



d 



2 r 



d(3 2 



D 
16i] 



-P 



We will use the WKB method to obtain a solution of the above equation [1 
Assuming that noise amplitude, D, is small, we seek the solution in the form 

P = exp[-S/D], S = W + DW X + D 2 W 2 + .... 



(20) 
2C 



(21) 



In the leading order in D one obtains a Hamiltonian-Jacobi equation for the function 
W: 



dW r] 
~~dz~ + 4 



dW 
drj 



1 2 



+ 



16r/ 



dW 
~df3 



0. 



The attached Freidlin-Wentzell Hamiltonian [20J is 

V 2 1 

4 P " ' 16r/' 



-A 



(22) 



(23) 
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where p v and pp represent the "momenta" canonically conjugated to the state variables, 
7] and (3 correspondingly: 

Pv = ~^, PP = W - (24) 
The solution of (|22|) is given by 



W 



, Mz') . ,,d{3{z') 



dz' - Ez, (25) 



where E = H is an integral of motion (see Refs. {TIJ|2n] for details). The Hamiltonian 
trajectories in the integrand of ({251) are f° un d from Hamilton's equations: 

d~z = 2 P ^ (26fl) 

£ = -z4 + W* (26c) 

^ = 0. (264 
dz 

These trajectories are subject to the boundary conditions: 

7,(0) =77o, ri{z)=V, P{0) = Po, P(z)=/3. (27) 
The solutions for the amplitude r](z) and the polarization angle (3{z) are given by: 

^ + 4r/ 2 (4 + z Pr?0 ) 2 
V{Z) = 64^ ' (28) 

ftz) = Po ~ arctan ( W) + arctan f 16 fa + '^ + 4 W ) . (29 ) 

\ P(30 J \ 877oP/30 

One can now find constants /Jgo and p,,o : 

,8^ yfm tan(/3 - | cos(/3 - /3 ) | 

P/30 — ± , 

z 

4 (-^±^|cos(^-A))|) 

Pt?0 = 



2 V^o 

From (J25j) - (j26dj) it follows that = £7z so that we have two solutions of Hamilton- 
Jacobi equation 



w± = 4 (r/ + r/ ±2 v ^ v ^|cos(/3-/3o)|) 



Statistics of a noise-driven Manakov soliton 



8 



The WKB approximation for the solution of (J2(J|) is written as a sum P(i],/3\z) 
P + (r),/3\z) + P-(r], /3\z), where P±(r), (3\z) is given by 



C 



det 



d 2 W± 
dqdq 



1/2 



exp 



Bpdz' 



x 



exp 



oqop 



exp 



(31) 



Here for the sake of compactness we introduced the following 2D vectors: q = (t],/3), 
<?o = (voiPo), P = (Pr/,Pi3)- Vector B represents the part of the advection vector which 
is proportional to D (with factor D omitted). In our case B = {3/4, cot(2(3)/8r]}. The 
integrals in the exponential are taken along the corresponding Hamiltonian trajectories 
d2HJ) ? (EH) , with the corresponding choice of sign in p v0 , pp Q . Factor C is added to provide 
the normalization of the PDF P(rj,(3\z). Each multiplier in (|31|) can be evaluated 
separately for both trajectories: 

d 2 W+\ 16 



det 



dqdq 



Tr d 2 H ^ _ 
dqdp I 



\f>« -Mi) 



and 



2 Z 

Bpdz' = J ^ + ^cot(2/3^ dz' = J Q + cot(2/3)/3^ dz' 







+ - In 


sin(2/?) 




Km) 


2 


sin(2/5 ) 



This yields 

P±{v,P\z)-- 



I sin(2/3) | 
sin(2/? )| 



exp 



L>2 



(77 + rjo ± | cos(/3 - (3 ) |) 



To obtain the normalization constant C we must integrate P + and P_ over 7] and /3 and 
insure that the sum of the two integrals evaluates to unity We perform the integration 
over p first. Since in the WKB approximation the noise intensity, D, is assumed to 
be small, the main contribution into the integral over f3 comes from the vicinity of the 
points where either cos(/? — (3$) or sin(/3 — (3q) vanish. This implies | sin(2/3)| ~ | sin(2/3 )|, 
and one obtains 

2tt 



P( V ,P\z)d(3 



16ttC rj 



exp 



Dz 



(V + Vo) 



Ir 



Dz 



(32) 



where we have used the identity / Q 27r cosh{x| cost\}dt = 27rl (x). The expression above 
is nothing else but a WKB approximation for the marginal PDF for the amplitude, 
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P(j]\z). Since the WKB approximation is valid under the assumption that D < 1 
(or, more precisely, Dz <C 1) one can use the asymptotic of the Bessel function for 
the large values of its argument and observe an excellent agreement with exact formula 
(I16|) obtained earlier. Another integration over rj allows one to obtain the normalization 
constant C. The final result for the marginal PDF P(t],/3\ z) reads: 



itDz t] V I sm(2/3 ) | 



cosh 



V^0| to o\\ 

1 cos(/3 — /3 Q )| 



Dz 



.(33) 



Note that this result is inapplicable for sin(2/3o) = 0. This is a direct consequence of 
the fact that the adiabatic perturbation theory fails for such values of polarization angle 
(see discussion at the end of Section EJ). 

4. Additive noise approximation 

Let us consider the case when the propagation distance z is small. We assume that the 
deviations of soliton parameters from their initial values are small as well. Then for 
generic set of Langevin equations (|A.1J) we may linearize the advection vector around 
the initial position q and assume that the perturbation functions g«(g, t) depends on 
the initial value of vector go rather than instantaneous value q (as in Appendix we omit 
vector notations for q and /). Under such assumptions Eq. (|A.l|) is transformed into a 
simple linear Langevin system with additive Gaussian noise: 

-fi(lo) + — [q-qo) + Si(z), (34) 

90 



with 



dz ' dq 

(si(z)sj(z')) = nij5(z - z'), (35) 

oo 

n i3 = 2GM = ~ Re J g*D gj dt. (36) 

— oo 

(The elements of the matrix G are derived in Appendix, see (jA.15|0 . The statistics of 
soliton parameters now become Gaussian. Similar approach is widely used to obtain 
approximate expressions for the variances of the soliton parameters in scalar noise-driven 
NLSE jlj. Without loss of generality one can assume that go corresponds to the steady 
state of the unperturbed system, i.e. fi(qo) = 0. The mean and variances of vector 
5q = q — go are now readily obtained in a tensor form: 

<Sq(z)>=0, (37) 

z 

< 5q{z) <g> Sq(z) >= 2 J dz exp[Az'}G (q ) exp[iV] , (38) 
df 



o 
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Now assuming that Qq = T = 0, and redefining phase a 
following expressions for the soliton jitters: 



< 5? >- 

< 5(p 2 > 



6 



288r7g 
D(Po,<Po), 



8r/ sin 2 (2/3 



D(n/2 + f3 , Vo ). i 



<5f3 2 >=—D(ir/2 + l3 ,(p ), 



a — 4:T]qZ we obtain the 

(40) 
(41) 

(42) 
(43) 

(44) 



and 



<6a 2 > = 2 F(/3q, <^ ) z + 2D 12 cot (2 fy) sin(2 <p ) z 2 rj 



128 



m + -2 



(45) 



Here we have introduced the function 

D(Po, Vo ) = D n cos 2 fa + D 22 sin 2 #, + D 12 cos(2p ) sin(2/3 ), (46) 

and D a p are the elements of the correlation matrix, see (j3J) • Formula 1)41)1 for the timing 
jitter was first obtained by Doktorov and Kuten [T£|. The expression for F(/3 ,ip ) in 
$ZE§ is rather involved: 



D(/? ,^ ) + — — {D u [129 + 4 7T 2 ] cos 2 (/3 )-36 (2 [D u + D 22 ] + 



48 7$ v ~ u ' ru/ ' 576t7o 
cos(2/3 )) + 9 cos(2/3 ) (l^ [ll + cot 2 (/3 )] + 4D 12 cos(2^ ) cot(2/3 )) + 
54 D n sin 2 (/? ) + 4D 22 [48 + vr 2 ] sin 3 (ft) + 9 tan 2 (/? ) sin 2 (/? ) + 
4 D 12 [12 + 7t 2 ] cos(2 ^ ) sin(2 (3 ) } . 

Eqs. ()40j) - ()42)) can be endowed with a simple interpretation. First one notices that a 
single soliton ansatz (0)),© can be transformed into a single soliton anzatz of the scalar 
NLSE via a unitary transformation: 



* = f{z) * 




cos/3(z) — sin/3(z) 
sin/3(^) cos/3(z) 







(47) 



(48) 



where Wo = 2rj sech[2r/(t — T)]exp[ia — 2 i(t — T)£] is a single soliton solution of the 
scalar NLSE. Since in the additive noise approximation we are interested only in 
small deviations of the soliton parameters from their initial positions we may assume 
T(z) ~ T(0). We can then observe that applying transformation T~ l leaves the l.h.s. of 
perturbed Manakov equation (JIJ) invariant. As for the noise in the r.h.s. it is transformed 
into n = T" 1 !!, which is again AWGN with the correlation matrix D = T~ l DT. 
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Therefore the variances of the amplitude, phase and frequency of Manakov soliton 
i n the additive noise approximation are equal to those of the scalar NLSE soliton 
but with the normalized noise intensity: D — > (-D)n = D((3 ,ip ). The variances of 
soliton jitters for scalar NLSE in this approximation are known jUEl- One can observe 
that formulae pn j) -(pi2 j) coincide with those given in [I] for amplitude, position and 
soliton frequency of a NLSE soliton provided that the noise intensity is now phase and 
polarization dependent D = D(^ ,ip ). 

This has a few interesting consequences. First we observe that for the isotropic 
noise when the correlation matrix D is proportional to the identity matrix, D = D I, 
amplitude, timing and frequency jitter do not depend on polarization or phase at all 
and we recover exactly the result for the scalar NLSE. In general case we must analyze 
the behavior of the function D{(3q, <fo). One can readily see that one of the local minima 
is given by: 



* a 1 
(p = - , p = - arctan 



2D 



12 



(49) 



£22 - D n 

The minimal value of the function D(/3q, (po) is 

Drmn = ~ | ~ - \ ^D\ 2 + (D 22 - D U )*. (50) 

Because of matrix D being positive definite det D > 0, and D m in is always non-negative, 
as it should be. If det D = 0, i.e. D n D 2 2 = Df 2 , then one can in principle dispose of 
the soliton jitters completely (in the additive noise approximation). In practice, one can 
make soliton jitters almost negligible, provided that the determinant of the correlation 
matrix D is close to zero. 



5. Conclusion 



To sum up, we have considered the non-Gaussian statistics of Manakov soliton and 
demonstrated the dependence of the variances of all soliton parameters on phase 
mismatch <po and polarization angle (3q. We derived the FPE describing the statistics 
of all soliton parameters. Using this approach we were able to work out the analytical 
expressions for the marginal PDFs for soliton amplitude 77, frequency £ and polarization 
angle (3. The statistics of the soliton frequency and amplitude appeared very much 
alike for scalar and vector NLSE in the case when vector AWGN is isotropic. This 
should not be considered as a surprise since for the isotropic perturbation of the 
Manakov system there exists a unitary transformation of the Manakov soliton into a 
NLSE soliton (localized in a single polarization) which leaves equation (jlj invariant [Hj. 
The situation however changes when one considers cross-correlations between noises in 
different polarizations. In the presence of a noisy perturbation with cross-correlations 
the symmetry related to the unitary transformation of a vector NLSE (space rotation 
and phase transformation) is broken. This means that though the relation between 
Manakov and NLSE solitons still exists the statistics of the soliton components becomes 
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polarization and phase dependent. Controlling the initial values of the polarization and 
the phase one can effectively manage the magnitude of each soliton jitter separately (in 
particular timing jitter). If the determinant of the noise correlation matrix is small one 
can effectively decrease the variance of a given soliton jitter and make it negligible. This 
is impossible for an isotropic noise and for a scalar NLSE. It means that the Manakov 
soliton is in general more robust against noisy perturbations than its scalar counterpart. 
The correlations between different polarization noise components can play a positive role 
in suppressing soliton jitters. 

The results obtained in the current paper rely on a series of assumptions. Firstly, 
the use of the adiabatic perturbation theory implies that we neglect the feedback of the 
linear radiation on the soliton. This is true for relatively small distances, where soliton 
parameters do not experience drastic deviations over the very short scales. Secondly, the 
perturbation theory developed here is unapplicable when the field polarization is close to 
linear (sin(2/3) ~ 0). Such cases should be considered separately (see |E] and references 
therein). And finally, we have used additive noise (or Gaussian) approximation to obtain 
the variances of the soliton jitters. This approximation works fairly well for relatively 
small distances (Dz <C 1) and gives reliable estimates for the variances. One should 
not be thinking, however, that the statistics of the tails of the PDF are Gaussian. In 
fact they are not as we have shown in Section EJ The FPE approach developed here 
may come handy when studying the statistics of the rare events forming the tails of the 
soliton PDF. 
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Appendix. Derivation of the Fokker-Planck equation 

The purpose of this Appendix is the derivation of the FPE which describes the 
PDF of the system with multi-component multiplicative noise. Suppose q s (z) = 
(ql(z) , . . . , q$f(z)) is a solution of set of N stochastic equations with M-component 
multiplicative white complex gaussian noise 



3C 




(A.l) 



— oo 



with 



n(t,z) 



(nti^z), . . .,n M (t,z) ) , 
(g iA (q,t),... } g iM (q,t)). 



(A.2) 
(A.3) 



and 



(n a (z,t)) = (n a (z,t)np(z',t')) = 0, 
(n a {z, t)n*Jz\ t')) = D a/3 5(z - z')5(t - t') 



(A.4) 
(A.5) 
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Here and in the sequel we will omit vector notations for the iV-dimensional vectors q 
and / in order not to confuse them with vectors gj and n since they have different 
dimensionality. 

We are interested in equation for the PDF of the solution of (jA.ljl . P(q,z). This 
PDF can be expressed in the form of 

P(q,z) = (6(q-q s (z))). (A.6) 

Differentiating the latter equation with respect to z and substituting the derivatives 
dq*(z)/dz one obtains 

oo 

BP 8 8 r 

= - £ dq-^ - £ R % J dt fe*(?> *) < % - 9'(z)Mt f z))}, (A.7) 

% ' — oo 

where the average in the latter term can be evaluated by virtue of Furutsu-Novikov 
formula [21] that in the case of vector white complex noise n(t, z) is generalized to 

<e[n,n>(M)>=D/§^\ , (A.8) 



5n*(t, z) 

where 0[n, n*] is an arbitrary functional of noise and matrix D has elements D a p (|A.5jl . 
Since 

5n* a {t,z) 



Sn}(t',z>) 
5n a (t, z) 



5 aP 5{z- z')5{t-t') , (A.9) 
0, (A.10) 



fin s (?) 1 1 

6(0)g k (q s (z),t), 9(0) = -, (A.ll) 



Sn*(t,z) 2 v ye> ^ w ' " vy 2 

(where in the latter formula we have used Eq. (|A.l|) and assumed symmetric Stratonovich 
regularization of the white noise) we arrive at 

(8(q - q s (z)) n(t, *)> = -i£ J-fi^ *)^(q» *)]■ (A.12) 

k 

Therefore for the PDF P(q, z) we arrive at the FPE equation: 

>) = - E |-l G <(«' p («- *)l + E aSs^M^. *)]■ < A - 13 > 

Here the components of the advection vector and the diffusion matrix (Gj and G^- 
respectively) are defined as 

oo 

G<(g) =Mq) + \ReJ2 I dt ||^ ( A - 14 ) 

fc -oo 

oo 

G lk (q) = hie j dtg*Dg k , (A.15) 

— oo 

where i, k = 1, . . . N. 
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